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In low-temperature pulsed growth two-dimensional islands form and coarsen into ~10 nm features. 
The islands produce well-defined displaced x-ray diffraction peaks due to relaxation of anisotropic 
surface stress of the (2x1) reconstruction with expansion and contraction present in orthogonal 
directions. We infer that the island distribution differs from continuous deposition, enhancing the 
population of size selected islands exhibiting anisotropic relaxation. The relaxation carries over into 
multilevel islands, suggesting that domains in subsequent layers form metastable stress domains. 



Surface stress plays a role in determining the equilib- 
rium structure of clean surfaces, as well as the evolution 
of structures during epitaxial growth. [I] Many surfaces 
exhibit anisotropic surface stress, which can influence the 
shapes of two dimensional islands, mound structures, and 
mesoscopic facets. [2HS] Even on surfaces where the low- 
est energy is a flat surface, crystal growth at low tem- 
perature can reveal anisotropic effects through the ap- 
pearance of low symmetry nonequilibrium features. Im- 
proved understanding of stress effects on surface growth 
may lead to approaches for use of elasticity as a tool for 
self-organization. 

On Ge(OOl) surfaces, neighboring surface atoms dimer- 
ize to form a (2x1) reconstruction to minimize the sur- 
face energy, but induce a large anisotropic stress. [Sj The 
stress is tensile parallel to the dimer bond, and com- 
pressive normal to it. Experimental support for stress 
anisotropy from the shape of islands, vacancy clusters 
and domain wall fluctuations have been reported. [H \5\ 
[3 H] In this work, we show that pulsed laser deposition 
(PLD) can produce nanostructures on Ge(OOl) surfaces 
which exhibit anisotropic lattice relaxation. Pulsed de- 
position appears to be advantageous for observation of 
these effects because the high instantaneous flux pro- 
motes the simultaneous nucleation of small islands that 
interact elastically at a very early stage and subsequently 
coarsen via strain-modified ripening. Observation of this 
process gives insight into the effects of anisotropic stress 
on the structures formed. 

The growth and annealing is monitored by real-time 
Grazing Incidence Diffraction (GID) around the (1 
0.05) s surface-sensitive reflection using synchrotron ra- 
diation with wavelength A=0.124 nm. Reflections are in- 
dexed in a surface unit cell with basis vectors related to 
the cubic axis as follows: [100] s = (l/2)[110] fcM/fc ; [010] a = 
(l/2)[110] 6u i fc ;[001], = [001]^. The sample prepara- 
tion and experimental setup details are given in a sup- 
plementary file. [5] 

The as-prepared vicinal Ge(001) surface consists of ~ 
50 nm wide terraces with the dimer orientation and stress 
tensor rotated 90° across each mono-atomic step, as il- 
lustrated in Fig. [ija). The large lattice relaxation con- 
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FIG. 1: GID evolution of (1 0.05) s reflection during sub- 
monolayer growth on Ge(OOl) by PLD at 100°C. (a) Schematic 
top view of small 2D dimer islands nucleated on a reconstructed 
surface containing degenerate (2x1) and (1x2) domains. Solid line 
represents the mono-atomic height step and big (red) circles are 
the dimer atoms. The arrows represent the elastic forces along the 
island periphery arising from the intrinsic surface stress anisotropy. 
The x-ray diffraction vector Q is oriented in the direction across 
the steps resulting from surface miscut. (b) The evolution of peaks 
E, B and C as a function of time. Dashed lines represent laser 
bursts: 1st, 150 pulses (0.3 ML), followed by 430 s recovery; 2nd, 
50 pulses to a total coverage of 0.4 ML followed by 1000 s recov- 
ery. The laser repetition rate is 10 Hz and the deposition rate is 
0.002 ML/pulse. (c) Diffraction profiles: (1) starting surface; (2) 
right after 1st laser burst; (3) before the 2nd laser burst; (4) after 
the 2nd laser burst and recovery. Note that the E peak is consis- 
tent with islands nucleated on Ta terraces while C corresponds to 
islands atop Tg terraces. 



centrates on small islands which are free to relax in two 
dimensions under the constraint of the covalent bonds to 
the substrate. The x-ray diffraction vector Q is oriented 
along the miscut, as shown. 

Fig. [ljb) shows the real-time GID intensity evolu- 
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tion near the (1 0.05) s reflection during low temper- 
ature growth at 100°C. Selected diffraction profiles are 
displayed in Fig. [l|c). Prior to growth, the diffraction 
profile shows only the main diffraction peak, labeled as 
B (Bulk), as expected for a well-ordered smooth surface. 
The sudden drop in intensity of peak B during the depo- 
sition burst is related to the formation of a high density 
of very small islands distributed over the surface. |10j [11] 

During the recovery time after the first deposition 
burst, a broad diffuse scattering background evolves into 
two displaced peaks E (Expanded) and C (Contracted) 
appearing on each side of peak B. Their integrated peak 
intensities without background subtraction are shown in 
Fig. [ijb) as a function of time. This relaxation be- 
havior is evidence for ripening as has been predicted for 
PLD.^TU] In the ripening process, small islands shrink 
until they disappear, while large islands grow at their 
expense. Similar effects have been observed for SrTiC>3 
homoepitaxy during PLD by specular x-ray scattering, 
which is sensitive to the island size and correlations. pTj 
Since we observe the effect with GID, our experiment is 
also sensitive to the lattice relaxation of the islands. 

The presence of two distinct diffuse peaks E and C 
suggests that the coarsened islands exhibit lattice relax- 
ation relative to the bulk value with both expansion and 
contraction. In particular, peak E is shifted down from 
the bulk position by 0.033 reciprocal lattice units (rlu) 
as indicated by curve 2, consistent with lattice expan- 
sion of islands atop T^ terraces in the direction parallel 
to Q. The relaxation is reduced slightly in subsequent 
curves and stabilizes at an average strain of about 2.5% 
expansion along the dimer rows. At 0.4 monolayer (ML) 
coverage, as shown in curve 4 of Fig. 1(c), the width of 
peak E along the radial direction (5h 0) s is 0.039 rlu. 
We interpret this width as domain broadening, so that 
it corresponds to an island size of L =10.3 nm along the 
dimer row direction. We conclude that the stabilization 
of both the island size and relaxation is due to the is- 
land size approaching a critical size, L c , where further 
ripening is inhibited due to an increase in strain energy. 
Thus, the largest island that can exhibit nearly complete 
relaxation is L c 10 nm. 

Peak C is weaker, presumably due to the fact that the 
density of dimer islands on the Tb terraces is less than 
on Ta terraces, resulting from the anisotropy in surface 
diffusion and bonding. [12] We also note that other effects 
can also influence the ratio of intensities, such as a small 
external stress on the surface, which can favor either the 
(2x1) or (1x2) at the expense of the other. [T3"] 

It is interesting to consider that all surfaces with 
anisotropic stress are in principle unstable to formation 
of elastic-stress domains. [U |2] For the case of Ge(001)- 
(2x1), Middel ct al. have found from measurements of 
the shape of vacancy islands that the stress anisotropy 
is (cr|| - a ± ) = 8.0 ± 3.0 eV/nm 2 . 5J Zandvliet et al. 
give step free energies at 100°C of F wa u— 0.048 eV/nm 
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FIG. 2: Evolution of the (1 0.05) s reflection in the multilayer 
growth regime, (a) The data is a continuation of the same de- 
position run shown in Fig. 1. The three curves correspond to 
deposited thicknesses of 0.4 (curve 4), 2.0 (curve 5) and 4.0 ML 
(curve 6), respectively. The curves are shifted vertically for clar- 
ity. Inset: schematic of a multilevel 2D dimer island illustrating 
the multilevel stress domain model. The wedding-cake-type stacks 
of (001) terraces are separated by mono-atomic height steps. The 
dimer rows and the resulting stress anisotropy rotates by 90° at 
each level, which is illustrated by the lines on each terrace and the 
arrows along the island periphery at each level, (b) Definitions of 
length L, width W, and island spacing d. as discussed in the main 
text. 

for Sa step edges and 0.104 eV/nm for Sb step edges. 
These values allow us to calculate the period A for both 
orientations of stress domains from 

x= Jjp_ ^(WC+i) (i) 



where C = (1 — v)(cr« — <7jJ 2 /27r/i, v is Poisson's ra- 
tio, and fi is Young's modulus. [U |^ Taking the domain 
fraction p to be 1/2 and the lxl lattice constant to be 
a = 0.40 nm, we have Xa = 22 nm and Xb = 86 nm. 
However, stress domains do not spontaneously form on 
the initial surface at the predicted length scales. We 
suggest that since there is a nucleation barrier to the 
formation of small islands the effect of the terrace steps 
keeps the adatom density on the terraces low enough to 
retard spontaneous nucleation of stress domains. In this 
case, stress domains should spontaneously form only if 
the terrace size can be made significantly larger, as has 
been observed for Boron doped Si(001)4l4 Below, we 
show from x-ray rocking scans that spatially correlated 
islands are formed in PLD, with spacings that closely 
agree with the values calculated from Eq. ([!]), lending 
credence to the idea that the pronounced strain relax- 
ation plays a role in stabilizing the structures observed. 
First, we will continue our discussion of the time-resolved 
data, now turning our attention to the > 1 ML regime. 

Fig. [2] is a continuation of the deposition run shown 
in Fig. [T] showing additional curves for growth thickness 
greater than 1.0 ML. Curve 4 is the same as in Fig. [ljc) 
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for deposited thickness of 0.4 ML, and curves 5 and 6 
are for thickness of 2.0 and 4.0 ML, respectively. The 
main trend in this growth regime is that peak B inten- 
sity continues to decrease and is eventually obscured by 
the diffuse intensity, while peaks E and C continue to in- 
crease. At 4.0 ML coverage (curve 6), diffuse peaks E and 
C suggest an average lattice relaxation of about 2.5% in 
expansion and 4.0% in contraction. After the completion 
of 4.0 ML growth, in-situ x-ray reflectivity shows a root- 
mean-square roughness of 0.28 nm (~2 ML), confirming 
that the growth is multilevel. 

This behavior is contrary to what would be expected 
if 2D islands grow beyond L c , since the lattice spacing 
on the interior of larger islands would be constrained to 
the bulk value. Moreover, in many cases, the in-plane 
lattice parameter is observed to oscillate because small 
islands become laterally constrained once they coalesce 
into a continuous layer. rTSHTT] . We do not observe os- 
cillatory relaxation related to 2D island coalescence be- 
cause multilevel growth takes place, i.e. nuclcation of 
the second monolayer occurs at ~ 0.5 ML and its stress 
anisotropy, which is rotated by 90° with respect to the 
layer beneath it, stabilizes the interior of the growing is- 
land. Each successive terrace in a given multilevel island 
is smaller than the terrace beneath it, thus a wedding- 
cake-type structure is formed, as illustrated in Fig. [2] 
The dimer direction and corresponding elastic force di- 
rection correspondingly rotate across each mono-atomic 
step gives rise to the alternating contributions to scat- 
tering intensity E and C. Thus, we infer the existence of 
metastable stress domains in multilevel growth. We re- 
fer to them as "metastable" since they are not believed 
to be lowest energy structures as compared to the two- 
level (2xl)/(lx2) stress domain envisioned by Alerhand 
et al. [2] This model most naturally explains the observa- 
tion that the strain relaxed peaks increase in intensity up 
to at least 4 ML deposited thickness without broadening 
significantly. 

Real space images of wedding-cake type multilevel is- 
lands have previously been observed in Scanning Tun- 
neling Microscope images after Ge growth on Ge(001) by 
Molecular Beam Epitaxy (MBE). [18] Our results are in 
striking contrast to the case of Ge MBE where layer-by- 
layer growth as a result of island coalescence precedes a 
transition to multilevel growth [91 [T5] As a result, for the 
case of MBE significantly thicker layers are required to 
observe the wedding-cake type structures. We explain 
these differences as follows: At a given growth tempera- 
ture the nucleation density in MBE is significantly lower 
than for the case of PLD because the peak incident par- 
ticle flux F is orders of magnitude lower, and the spacing 
between nuclei varies as (D/F) 1 / 6 , where D is the surface 
diffusion constant. [19] Island sizes may exceed L c , and 
the monomer density on the surface is subsequently very 
low, so that nucleation of new small islands is kinetically 
blocked. In contrast, for PLD, due to the high island 



nucleation density, neighboring islands start to interact 
with each other through elastic interaction at an early 
stage, thus favoring size selection and ordering. Viewed 
another way, we can also say that the islands are at the 
right length scale to break the nucleation barrier for the 
formation of stress domains. Moreover, when strain fields 
of neighboring islands overlap, coalescence can be inhib- 
ited in favor of multilevel growth. This effect has been 
observed in Monte Carlo simulations. [20] Alternate mod- 
els based on standard Ostwald ripening do not predict the 
formation of regular arrays of equally spaced islands. [2"U] 
The large lattice relaxation during Ge(001) homoepi- 
taxy illustrated in Figs, [I] and [2] is consistent with a 
recent observation of 3.6% contraction in the dimer bond 
direction near the Sa step edge on Si(001) by non-contact 
atomic force microscopy at 5 K.[3T] It is found that the 
elastic relaxation extends to about 4.0 nm away from 
the step. If we consider a 2D island edge instead of a 
straight step edge, the corresponding critical diameter 
above which the islands can no longer be fully relaxed 
would be around L c — 8.0 nm. This predicted island size 
agrees well with our observation of ss 10 nm island size 
at 0.4 ML. 
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FIG. 3: RHEED patterns obtained along the [100] s azimuth of (a) 
as-prepared vicinal Ge(001) surface; (b) 4.0 ML thick layer grown 
at 100°C; (c) 5.2 ML thick layer grown at 250°C. The electron 
energy used is 20 keV with incident angle of 2.2°, 3.0° and 2.8°, 
respectively. 



The contraction resulting in the C peak in our data 
is associated with the dimer bond formation, and the 
expansion in the orthogonal direction resulting in the E 
peak is driven by tilting of the dimer, since the lower 
dimer atom tends to spread out and push the atoms be- 
neath it laterally. [5] We calculated subsurface displace- 
ments by minimizing elastic energy with the Keating 
model. Displacements calculated for a perfect 2x1 sur- 
face structure show elastic relaxation to at least six layer 
beneath the surface, consistent with previous experimen- 
tal and theoretical results. [55] Calculations for a mini- 
mal four dimer island suggest that the dimer layer and 
at least two layers beneath the dimer island have large 
displacements (> 1%), and should contribute to both 
E and C peaks. The calculation also shows that large 
displacements propagate into the surrounding substrate 
surface layer. This is consistent with the idea that strain 
relaxation may influence attatchmcnt/detachment pro- 
cesses at the island edges, as well as introducing a bias 
to diffusion near the islands. These are possible atomic 
scale mechanisms to explain the observation that island 
sizes saturate at L c during ripening. Both processes have 
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also been suggested to play a role in the ordering of self- 
organized quantum dots. [23ll24] 

To link the film morphology with results found in the 
literature, Reflection High Energy Electron Diffraction 
(RHEED) patterns were recorded along the [100] s az- 
imuth before and after growth. Fig. [3] shows several 
examples. On the clean surface before deposition, in 
(a) bright spots coexist with half-order reconstruction 
streaks, indicating a smooth and well-ordered Ge(001)- 
(2x1) reconstructed surface. 25J After finishing growth, 
in RHEED patterns (b) and (c) streaks broaden, half- 
order streak intensities decrease while diffuse scattering 
is observed to increase, consistent with the surface rough- 
ness deduced from x-ray reflectivity data. In addition, 
intensity modulations become visible along the integer 
diffraction streaks. This behavior has previously been 
correlated with the formation of multilevel 2D islands 
on the growing surface. [26l [27] i.e., the growth mode 
in which several levels separated by mono-atomic height 
steps are exposed on the surface. 
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FIG. 4: (a)-(b) Reciprocal space maps near (1 0.05) s of the 
sample deposited at 100°C: (a) prior to growth and (b) after 4.0 
ML growth. Plots (c) and (d) show results for a new sample with 
a deposition temperature of 250° C and subsequent annealing at 
650° C for 400 s and cooling down to room temperature. Rocking 
scans have been converted from angles to reciprocal lattice units. 
The line profile in (d) is along the transverse direction (0 8k 0) s for 
peak E, as indicated by the vertical white line in (c). The dotted 
lines are fitted Gaussian profiles. 

In order to obtain additional information about the 
multilevel island structure, we performed rocking scans 
on the sample about its surface normal. A linear position 
sensitive detector was used so that 2D reciprocal space 
maps can be rapidly obtained with a single scan. Wc 
first describe how lattice relaxation, island size broad- 



ening and island correlation affect the GID data: (i) In 
reciprocal space near (100) s , radial (Sh 0) s intensity 
profiles are sensitive to the strain relaxation parallel to 
the surface. The peak position indicates the average re- 
laxation, while the peak width can be limited by the dis- 
tribution of lattice constants or by domain size broaden- 
ing. In the latter case, the width is approximately a/L, 
where the surface unit cell constant is a — 0.400 nm and 
L is the island size along Q as shown in Fig. [2j (ii) 
Transversely, (0 5k 0) s diffuse scattering gives informa- 
tion about the extent of correlations laterally across the 
surface. In cases where islands are randomly distributed 
on the surface, the diffuse scattering will be centered at 
k = and its width will be approximately a/W, where 
W is the average island size, (iii) In cases where island 
positions are spatially correlated, satellite peaks occur 
near Ak = ±a/d, where d is the average island-island 
spacing. [28] 

Fig. |I](a)-(b) shows detailed reciprocal space maps at 
the (1 0.05) s position of the sample (a) prior to growth, 
(b) after 4.0 ML growth at 100°C. In Fig. ga), a high 
intensity peak B indicates a well-ordered starting sur- 
face. In Fig. |4|b), peak B is no longer visible, due to 
the surface roughening. Two satellite peaks are observed 
vertically on each side of peak E, near Ak = —0.018 and 
+0.021 rlu. This indicates a strong spatial correlation 
of multilevel 2D islands along the dimcr bond direction 
with a uniform island spacing of d ~ 20 nm. The effect is 
anisotropic because no correlation is observed in the or- 
thogonal direction. This value of d is in close agreement 
with the length scale calculated for stress domains, 
and is reasonable if we assume that the length scale was 
set during the submonolayer stage of growth, given that 
the islands are separated by Sa steps, as shown in Fig. 
[2jb). We also find that L ss 16.9 nm. To complete our 
description of Fig. |4|b) , we note that the widths of peak 
C indicate an average island size of L rs 15.1 nm along 
the dimer direction, and W ~ 12.3 nm in the transverse 
direction. Since peak C only increases rapidly for depo- 
sitions > 1 ML, we interpret these widths as being pre- 
dominantly relevant to domains within multilevel islands 
atop Tj\ terraces. 

Additional information about surface kinetics can be 
obtained by annealing the metastable multilevel struc- 
tures discussed so far. Fig. [4] (c) shows an additional re- 
ciprocal space map on another sample deposited at 250° C 
with a total thickness of 5.2 ML, followed by annealing 
at 650° C. It is noted that the new sample shows almost 
the same peak E and C positions and widths as the one 
deposited at 100°C. This shows that the island size again 
saturates near L c rs 10 nm for submoloayer deposition, 
and that it is a true saturation rather than simply a re- 
sult of the kinetics slowing down as the island size in- 
creases because the coarsening process is much faster at 
250°C. Upon annealing peak positions stay unchanged. 
However, the profile of peak C becomes more compact 
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along the k direction, while peak E becomes more com- 
pact along the radial h direction. Specifically, peak E 
changes width along h, corresponding to a change in the 
island size along the dimer row direction from L rj 12.5 
to 20.0 nm. Peak C also becomes more compact in k, 
suggesting that W changes from < 10 to 16.7 nm. The 
satellite peaks for diffuse scattering E shown in Fig. |4|d) 
after annealing are located at Ak — ±0.019 rlu, showing 
that the island spacing along the dimer bond direction 
has changed slightly from D 18.8 nm to 21.0 nm, in- 
dicating that minimal ripening has occurred and hence 
the majority of the mass transport is related to transport 
within the multilevel islands. 

We have annealed stepwise to higher temperatures un- 
til the metastable islands relax away and the surface re- 
turns to a nearly flat state at 800° C. In a separate set 
of experiments, we also confirmed that the surface phase 
transition occurs at sa 780° C. The transition is charac- 
terized by reversible broadening of the (3/2 0.05) re- 
flection, which disappears entirely above the transition 
temperature, and loss of the ordered terrace structure, 
as confirmed by monitoring the peak splitting in k scans 
at (0 1 0.05). Our observations are consistent with the 
suggestion of Zandvliet et ah, [25] that the transition is 
due to a roughening transition where the step edge free 
energy F wa u is reduced to zero rather than by breaking of 
dimer bonds. In this case, Eq. [T] predicts a length scale 
of A a: 7 nm with no anisotropy, indicating that small 
3-4 nm islands may still be present on the surface during 
annealing at 800° C, but these are not directly detected 
in our measurements. 

In conclusion, we have observed two main effects: (i) 
anisotropic relaxation in submonolayer deposition dur- 
ing pulsed growth with saturation of island sizes at ^10 
nm, and (ii) persistence of anisotropic relaxation during 
multilevel growth and subsequent annealing to moder- 
ate temperatures. The results show that high density 
nucleation followed by evolution through a ripening pro- 
cess leads to a narrow distribution of size-selected islands, 
which can also be viewed as breaking the kinetic barrier 
to the formation of stress domains. 

The authors acknowledge Christie Nelson and Steve 
Lamarra for experimental assistance with the work done 
at the NSLS X21 beamline. Research supported by the 
U.S. Department of Energy, Office of Basic Energy Sci- 
ences, Division of Materials Sciences and Engineering un- 
der Award No. DE-FG02-07ER46380. Use of the Na- 
tional Synchrotron Light Source was supported by the 
U.S. Department of Energy. Development of the capa- 
bility for in-situ x-ray analysis of PLD was supported 
by the National Science Foundation under DMR-0216704 
and DMR-0348354. 



[1 

[2 

[3 

[4: 

[5 

[6 
[7 

[8 

[9 
[10 

[11 

[12 
[13 
[14 

[15 
[16 

[17 

[18 
[19 



Present address: Nanjing University of Aeronautics and 

Astronautics, Nanjing, China. 

rheadrick@uvm.edu 

Ibach, H. The role of surface stress in reconstruction, epi- 
taxial growth and stabilization of mescoscopic structures. 
Surf. Sci. Rep. 29, 193 (1997). 

Alerhand, O. L., Vanderbilt, D., Meade, R. D. & 
Joannopoulos, J. D. Spontaneous formation of stress 
domains on crystal surfaces. Phys. Rev. Lett. 61, 1973 
(1988). 

Kochanski, G. P. Step-step interactions due to 
anisotropic surface stress. Phys. Rev. B 41, 12334 (1990). 
Li, A., Liu, F. & Lagally, M. G. Equilibrium shape of 
two-dimensional islands under stress. Phys. Rev. Lett. 
85, 1922 (2000). 

Middel, M. T., Zandvliet, H. J. W. & Poelsema, B. Sur- 
face stress anisotropy of Ge(001). Phys. Rev. Lett. 88, 
196105 (2002). 

Zandvliet, H. J. W. The Ge(001) surface. Phys. Rep. 
388, 1 (2003). 

Wu, F. & Lagally, M. G. Ge-induced reversal of surface 
stress anisotropy on Si(001). Phys. Rev. Lett. 75, 2534 
(1995). 

Zandvliet, H. J. W. & Poelsema, B. Determination of 
surface stress anisotropy from domain wall fluctuations. 
Phys. Rev. B 59, 7289 (1999). 

See EPAPS Document No. [to be inserted by publisher] 
for sample preparation details and additional figures. 
Vasco, E., Polop, C. & Sacedon, J. L. Preventing 
kinetic roughening in physical vapor-phase-deposition 
films. Phys. Rev. Lett. 100, 016102 (2008). 
Ferguson, J. D., Arikan, G, Dale, D. S., Woll, A. R. 
& Brock, J. D. Measurements of surface diffusivity and 
coarsening during pulsed laser deposition. Phys. Rev. 
Lett. 103, 256103 (2009). 

Voigtlander, B. Fundamental processes in Si/Si and 
Ge/Si epitaxy studied by scanning tunneling microscopy 
during growth. Surf. Set. Rep. 43, 127 (2001). 
Men, F. K., Packard, W. E. & Webb, M. B. Si(100) 
surface under an externally applied stress. Phys. Rev. 
Lett. 61, 2469 (1988). 

Ermanoski, I., Bartelt, N. C. & Kellogg, G. L. Self- 
assembly of defect-free nanostripe arrays on B-doped 
Si(001). Phys. Rev. B 83, 205432 (2011). 
Massies, J. & Grandjean, N. Oscillation of the lattice 
relaxation in layer-by-layer epitaxial growth of highly 
strained materials. Phys. Rev. Lett. 71, 1411 (1993). 
Eymery, J., Daudin, B., Cunff, D. B., Boudet, N. & 
Tatarenko, S. Anisotropic relaxation during the first- 
stage of the growth of ZnTe/(001)CdTe strained layers 
studied by reflection high energy electron diffraction. J. 
Appl. Phys. 66, 3456 (1994). 

Hartmann, J. M., Arnoult, A., Carbonell, L., Etgens, 
V. H. & Tatarenko, S. Reflection high-energy electron 
diffraction measurement of lattice-parameter oscillations 
during homoepitaxial growth of CdTe. Phys. Rev. B 57, 
15372 (1998). 

Van Nostrand, J. E., Chey, S. J. & Cahill, D. G. Low 
temperature growth morphology of singular and vicinal 
ge(001). Phys. Rev. B 57, 12536 (1998). 
Pimpinelli, A. & Villain, J. Physics of Crystal Growth 



6 



(Cambridge University Press, 1998). 

[20] Scholl, E. & Bose, S. Kinetic monte carlo simulation 
of the nucleation stage of the self-organized growth of 
quantum dots. Solid-State Electron. 42, 1587 (1998). 

[21] Naitoh, Y., Li, Y. J., Nomura, H., Kageshima, M. & Sug- 
awara, Y. Effect of surface stress around the Sa step of 
Si(001) on the dimer structure determined by noncontact 
atomic force microscopy at 5 K. J. Phys. Soc. Jpn 79, 
013601 (2010). 

[22] Pedersen, J. S. Surface relaxation by the keating model: 
a comparison with ab-inito calculations and x-ray diffrac- 
tion experiments. Surf. Sci. 210, 238 (1989). 

[23] Pan, E., Zhu, R. & Chung, P. W. On the correlation 
between the self-organized island pattern and substrate 
elastic anisotropy. J. Appl. Phys. 100, 013527 (2010). 

[24] Meixner, M. & Scholl, E. Kinetically enhanced correla- 
tion and anticorrelation effects in self-organized quantum 



dot stacks. Phys. Rev. B 67, 121202(R) (2003). 

[25] Horn, K. M., Chason, E., Tsao, J. Y., Floro, J. A. & 
Picraux, S. T. Oxygen roughening of Ge(001) surfaces. 
Surf. Set. 320, 174 (1994). 

[26] Xue, C, Xiao, H. Z., Hasan, M. A., Greene, J. E. 
&: Birnbaum, H. K. Critical epitaxial thickness for 
low-temperature (20-100°C) Ge(100) 2x1 growth by 
Molecular-beam Epitaxy. J Appl. Phys. 74, 2512 (1993). 

[27] Bratland, K. A. et al. Mechanism for epitaxial break- 
down during low-temperature Ge(001) molecular beam 
epitaxy. Phys. Rev. B 67, 125322 (2003). 

[28] Guinier, A. & Fournet, G. Small-angle scattering of x 
rays (Wiley, 1955). 

[29] Zandvliet, H. J. W. Determination of Ge(001) step free 
energies. Phys. Rev. B 61, 9972 (2000). 



